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1. Introduction 

Recent experiments on quantum gases in one-dimension (ID) [H, El El IH [5] provide an 
exciting way to test ID many-body physics extensively studied in the literature via the 
Bethe ansatz (BA) [6], effective field theory [7] and other theoretical methods. Novel ID 
many-body physics such as the phenomena of spin-charge separation (HI |9j HU] , universal 
Luttinger liquid thermodynamics [HI [121 1131 EEH US] and quantum criticality, are quite 
different from higher dimensional physics. ID many-body systems also provide insights 
into higher dimensional physics. For example, theoretical predictions for the existence of 
a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) [16] like pairing state in the ID interacting 
Fermi gas have emerged [131 HZl HHJ UHl [20l [2U [22l [231 12H I25l [261 12Z] . 

A scheme for mapping out the physical properties of homogeneous systems by the 
inhomogeneity of the trap has been successfully applied to ultracold atom experiments 
[5], since the thermodynamics of interacting fermionic systems can be measured 
explicitly and the interaction between fermions can be tuned precisely via a broad 
Feshbach resonance. In particular, fermionic alkaline-earth atoms display an exact 
SU(N) spin symmetry with N = 21 + 1 where I is the nuclear spin [28]. For 
example, a recent experiment dramatically realized the model of fermionic atoms with 
SU (2) x SU(6) symmetry where electron spin decouples from its nuclear spin 1 = 5/2 for 
171 Yb [29] . Such fermionic systems with enlarged SU (N) spin symmetry are expected to 
display a remarkable diversity of new quantum phases and quantum critical phenomena 
due to the existence of multiple charge bound states [301 El|- The strongly attractive 
mult i- component Fermi gases with higher spin symmetries have also been investigated 
via the Bethe ansatz method [321 EH El]. These developments provide an exciting 
opportunity to explore the universal thermodynamics and quantum critical behavior of 
strongly interacting fermions with high spin symmetries in ID. 

The Bethe ansatz has proven to be a powerful method to study ID quantum gases. 
This method was applied to the ID Bose gas by Lieb and Liniger [35J and to the spin- 1/2 
Fermi gas by Yang [36J and Gaudin [37]. After that, the multi-component Fermi gas was 
studied by Sutherland [38J. The study of the thermodynamics of the attractive Fermi 
gas was initiated by Yang [39] and Takahashi [40J. Moreover, the thermodynamics 
of ID Fermi gases is also widely studied for its rich physics. For spin-1/2 fermions 
with attractive interaction, the phase diagrams and Luttinger liquid physics have been 
investigated by using both analytical [131 E] an d numerical methods [TBI [T7]. The 
key features of this T = phase diagram were experimentally confirmed using finite 
temperature density profiles of trapped fermionic 6 Li atoms [5]. 

For the strongly attractive spin-1/2 Fermi gas at finite temperatures, the 
thermodynamics of the homogeneous system is described by two coupled Fermi gases 
of bound pairs and excess fermions in the charge sector and ferromagnetic spin-spin 
interaction in the spin sector [131 IT4] . Spin fluctuations are suppressed by a strong 
effective magnetic field at low temperatures. However, for the repulsive case, the 
situation is quite different since spin fluctuations play a dominant role in the low 
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energy physics. Spin fluctuations are described by the SU{2) antiferromagnetic spin- 
spin Heisenberg system and thus the coupling of charge and spin parts make the 
thermodynamic Bethe ansatz (TBA) equations more complex. So far we know that the 
Wiener-Hopf method can be applied here to deal with the spin part of the TBA equations 
for weak magnetic field. This method was used by Mezincescu et al. to investigate the 
thermodynamics of Heisenberg spin chains with SU(2) and SU(3) symmetries [HJ @2]. 
Recently, Lee et al. [15] applied the Wiener-Hopf method to spin-1/2 repulsive fermions 
with SU(2) symmetry and studied the low temperature thermodynamics and correlation 
functions via the Bethe ansatz method analytically. The Wiener-Hopf technique was 
applied to get the solutions of the TBA equations in the spin sector. Then Sommerfeld's 
expansion was utilized in the low temperature limit with strong interactions to get the 
multi-component Tomonaga-Luttinger liquid (TLL) form of the free energy, which shows 
a spin-charge separation with the central charge of the spin and charge parts both equal 
to 1. 

Spin-charge separation is a hallmark of ID many-body physics. It is a universal 
feature that interacting particles "split" into spins and charges as temperature tends 
to absolute zero temperature. The collective excitations with only spin or charge are 
called spinon and chargon/holon (the antiparticle of chargon), which have different 
velocities. This behavior was studied in different kinds of materials [HI El HQ] arid recently 
observed in experiments, for example, of ID metallic wires on surfaces, ID organic wires, 
carbon-nanotubes, quantum wires in semiconductors, and other types of ID systems. 
ID quantum systems with high spin symmetries will give diverse magnetism and exotic 
properties for these kinds of gases. However, such systems are more difficult to solve 
[43l 1441 145] due to the complicated magnetic ordering. The three-component Fermi gas 
has U(l) x SU(3) symmetry which leads to two sets of spin waves. It is very interesting 
to see how the low temperature thermodynamics of such a gas naturally separates into 
free Gaussian field theories for the U(l) charge degree of freedom and the two spin 
sectors. 

In this paper, ID fermions with U(l) x SU(3) symmetry is studied by using 
the Wiener-Hopf method. The result shows that the system with strong repulsive 
interactions at low temperatures has the central charge C c = 1 for the charge sector 
and central charge C s = 2 for the two spin sectors. We also derive the leading order 
finite temperature corrections of the free energy which is consistent with expectations 
from conformal field theory (TTJ [12j 06j H7] . In addition, the polylog function is used to 
derive the equation of state in a wider temperature regime. 

This paper is organized as follows. In Section II, the BA and TBA equations 
are presented. In Section III, the Wiener-Hopf method is applied to solve the TBA 
equations for the spin sectors. In Section IV, the thermodynamics of the system is 
derived via Sommerfeld's expansion for low temperatures and strong interactions. We 
also show that the expression for the free energy at low temperatures is in agreement 
with conformal field theory and derive the spin and charge velocities. Section V is 
reserved for the conclusion and discussion. Some details are given in the appendices. 
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2. The Model and thermodynamic Bethe ansatz 

Consider a ID system of N fermions with mass m and a spin independent 5-function 
potential. The fermions can occupy three possible hyperfme levels (|1), |2), and |3)) 
with particle numbers N , N 2 , and iV 3 , respectively. They are also constrained to a 
line of length L with periodic boundary conditions. This system is described by the 
many-body Hamiltonian J38J HQ] 



h 2 N d 2 

^0=-r-Ez2+fflfl E S( Xi -Xj)+E g . (1) 
£Ul i=1 CXj l<i<j<N 

Here E z = E*=i N% is the Zeeman energy. The spin independent contact interaction 
giD exists between fermions with different hyperfme states so that the number of 
fermions in each spin state is conserved. It is positive for repulsive interaction and 
negative for attractive interaction. For simplicity, we set h = 2m = 1. They can be 
reintroduced when necessary. It is possible to tune scattering lengths between atoms in 
different sublevels to form nearly SU(3) degenerate Fermi gases via a broad Feshbach 
resonance. 

The Hamiltonian in equation (TjQ) exhibits a symmetry of U(l) x SU(3), where U(l) 
and SU(3) are the symmetries for charge and spin degrees of freedom, respectively. 
As mentioned earlier, this model was solved via a nested BA [381 HO]- The energy 
eigenspectrum is given by E = Y.f=i k 2 in terms of the quasimomenta {kj} satisfying 
the BA equations SO] 



e 



e= \ kj - A e - ic/2' 

SL Aj - kj + ic/2 _ ^ A e -A a + ic ^ A e - \ rn - ic/2 
^ A, - k 3 - ic/2 " M A e -A a - ic ii x A e -X m + ic/2 ' 

J| A m -A, + ic/2 = _ Q X m -Xfs + ic (2) 
£ =1 X m — An — ic/2 „ =1 X m — Xp — ic 

Here j = 1, . . . , N, £ = 1, . . . , Mi and m = 1, . . . , M2 with quantum numbers 
Mi = N 2 + N 3 , and M 2 = N 3 . We define the interaction strength as c = mgio/h 2 > 
since we only consider the repulsive case. The parameters {kj} are quasimomenta 
and {Ai, X m } are the rapidities that characterize the internal hyperfme spin degrees of 
freedom. In the thermodynamic limit, N, L — > 00 with the ratio for linear particle 
density n = N/L kept finite. Notice that there are two classes of strings. Both sets of 
solutions can be expressed in terms of the string hypothesis 

A"* a = A"-I(n + l-2a)i|c|, a = l,...,n (3) 
2 

Xf =A?-~(n + l-2/3)i|c|, /3 = l,...,n (4) 

where n is the length of each string, j labels each individual string, and A™ and A™ 
are the real parts of the A and A strings, respectively. To distinguish the rapidities in 
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different spin levels, define A( r ) as 

A:=A«, A:=A< 2 >. (5) 

For repulsive interaction, the quasimomenta {kj} are real, while {A^} form 
complex spin strings which characterize the spin wave fluctuations at finite temperatures. 
Without loss of generality, we consider the case where equal Zeeman splitting occurs 
between fermions of different spin to simplify our analysis, i.e., Hi — H 2 — H. The 
equilibrium states are then determined by minimizing the Gibbs free energy, which gives 
rise to a set of coupled nonlinear integral equations - the TBA equations, which are 

00 m 

e(k) =k 2 -fi-H-Tj2an* hi(l + e~*» (fe)/T ), 

n=l 

^(AW) =nH- Ta n * ln(l + e~^ T ) + T]TT mn * ln(l + e -^^ T ) 

m 

m 

^>(A< 2 >) =nH + T^2 T mn * ln(l + e~^^ T ) - T £ S mn * ln(l + e'^^/ft) 

m m 

Here e(k) denotes the dressed energy and <f>^ (r = 1,2) are associated with densities 
of strings with length-n in \^ parameter spaces. * represents the convolution integral, 
i.e. (a * b)(x) = J a(x — y)b(y)dy, and the function a m {x) is given by 
/ \ _ f mc 

am[X) ~2^(mc/2y + X 2- [V 
The functions T mn and S mn are defined as 

rp / x = f a m+n (x) + 2a m+n - 2 (x) + • • • + 2a\ m - n \ +2 {x) + a\ m - n \(x), for m ^ n\ . . 
\ 2a 2 (x) + 2a 4 (x) + . . . + 2a 2n _ 2 (x) + a 2n (x), for m — n. 

g ( x ) = I a m+n-i(x) + a m+n _ 3 (x) + . . . + a| m _ n | +3 (x) + a\ m - n \ +1 (x), for m ^ n- 
\ ai(x) + a 3 (x) + . . . + a 2n - 3 (x) + a 2n -i(x), for m = A. 

An alternative set of TBA equations are 

e(k) = k 2 - fi-T Sl * ai * ln(l + e - £{k)/T ) 

- Ts, * ln(l + e^^ T ) - Ts 2 * ln(l + e ^^ T ), 
J>?(\W) = Ts * In (l + e ^ A(1) )/ T ) - Ts * In (l + e -< xW » T ) 

-Ts*\n(l + e-^^^ 



with the limit 



^)(AW) = Ts * In (l + e<^ A(1) )) + Ts * In (l + e O* (1) )) , 
#(A (2) ) = Ts * In (l + e ^( A(2) )/ r ) - Ts * In (l + 
0(f) ( A ( 2 )) = Ts * In + e ^ (2) )) + Ts * In (l + e O* (2) )) . (10) 

li m ^ lA ; = # (r = 1, 2) (11) 
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and 

= 2 C COsh (vrx/c)' (12) 

1 cosh (ttx/3c) 

~ V3c cosh (tts/c) ' ^ j 

, , 1 sinh (irx/3c) . , s 

s 2 (x) = —= — , - \ ' / . 14 
v ' y/3c smh (ttx/c) 

The two sets of TBA equations are interchangeable via Fourier transformation. 
3. Wiener-Hopf Solution For The Spin Part 

The string part is the most difficult part of the TBA equations to solve. It consists of an 
infinite number of string functions 0^(A) which are related through the set of equations 
(JO]) or (ITU]) . These coupled nonlinear integral equations cannot be solved analytically 
in most cases. However, in some special cases the equations can be simplified. When 
T is very low (T 1), these equations reduce to a set of linearly coupled equations. 
Moreover, when c is very large (c> 1), Taylor expansion can be applied to solve this 
set of equations. In the following, we will consider the case where T 1 and c ^> 1. 
Observe from the second set of TBA equations ([TO]) that > for n > 1 at low 

temperatures. The functions s(A) and In ^1 + e^™' > ( A )/ T ^ are always greater than zero 
in their entire domain, so the convolution is always greater than zero. This positivity 
condition implies that the function Tin ^1 + e - ^™'^/ 7 ^ — > for T — > 0. Thus all the 

spin string functions could be neglected except the n = 1 function ^^(A) in the first set 
of TBA equations. For T < 1, only the lowest strings are therefore left in the strong 
coupling limit c> 1, and we can rewrite the term as 

Ta x * In (l + e~ £ik)/T ) « 27rPa 1 {k), (15) 



where 

The first set of TBA equations ([6]) can thus be simplified as 



/oo , , 

dkln(l + e- £{k)/T ) . (16) 



e (k) =k 2 -n-H-T ai * ln(l + e"^ 1)(fc)/T ) 



01 (A) = H- 2nPa 1 (X) + Ta 2 * In ( 1 + e - ^ {A)/T J - Ta x * In M + e - ^ (; 

0f)(A) = H + Ta 2 * ln(l + e -*£ 2) W/ r ) - Ta x * ln(l + e^^). 

where we have switched to a common spin variable A to represent both spin spaces A^ 
and \( 2 \ Taking the limit T — > yields 

e (k) =k 2 -^-H + a 1 * (k) , (20) 
<tf } (A) = H - 2uPa l {k) - a 2 * ^ (\) + o x * </>i 2) ~(A), (21) 
0( 2 ) (A) = H - a 2 * 0i 2) " (A) + oi * <rf 1)_ (A). (22) 
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We have decomposed the functions as 

^ r) (A) = # + (A) + 4 r) -(A), (23) 

where 

fc M+m_ / ^W, for0S r) (A)>O; 

0, for 0f' } (A) < 0. 



W +(A) = ^ ^ ^, - ^ - u, (24) 



^(A) = Pi) for t (A)> ° ; (25) 
' ^ <rf r) (A), for <rf r) (A) < 0. 



After taking the Fourier transfers and going through some manipulation, equations 
( l2~TT) and ( 1221) can be written as 

<# } (A) = ii - 2ttP Si (A) + * ^ 1)+ (A) + g * # + (A), 

0S 2) (A) = H- 2ttPs 2 (A) + /i * <M 2)+ (A) + 5 * 0S 1)+ (A), (26) 

where 

h(X) = s 1 *a 1 (X)-s 2 (X), (27) 
ff (A) =s 2 *a 1 (A)-s 1 (A). (28) 
Both equations in ([2"B"j) are the T — )■ limit cases of 

0«(A) = if - 2ttP Si (A) + h * Tln(l + e ^ 1)(A)/T ) + <? * Tln(l + e<^ )(A)/T ), 
0f)(A) = if - 2ttPs 2 (A) + h * Tln(l + e^ 2)(A)/T ) + * Tln(l + e^ 1)(A) (2» 

In order to find a relationship between the dressed energies and the temperature, 
we write <j)f\X) as 

#(A) = 0M(A)+r^(A), (30) 

where the first term 0( r )(A) is the term of zeroth order when T = and second term 
(A) is the first order correction to the limit T — > 0. Substituting this equation into 
(T2T)|) and using the expressions in f[2"B"j) for (fi^(X), the equation for r/ r )(A) becomes 

^D(A) = * [Tin (l + e^ (1) ( A ) + ^ (1) ( A ))/ T ) - 0«+(A)" 

+ 0* [Tln(l + e^ (2) ( A ) + " (2) ( A ^ T ) -</> (2)+ (A)] , (31) 

77^ (A) = h* [Tin (l + e ^ (2) ( A ) + " (2) ( A »/ T ) - 0< 2)+ (A)] 

+ 0* [rin(l + e^ (1) W + " (1) ( A »/ T ) -0 (1)+ (A)] . (32) 

Since <// r )(A) is an even function and A = Aq gi ye the zero points, the above 
equations for 77^ (A) can be simplified to become 

J\\\>\ 

+ f , , 2) MA - X')^(X')dX' + e£\X) + E® (A), (33) 

./|A|>A 

^ (2) (A)^ i _^h(X-X')rj^(X')dX' 

(i) h(X - X') V ^(X')dX' + 4 2) (A) + EM(X), (34) 



|A|>A< 2) 



|A|>A« 
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where the approximation 

/* [Tln(l W^ M(A)+??(r){A))/T " 



|A|>A& 



(r) 



(1)+ (A) 

dX'f(X-X') [Tln(lW^ )(A)+ ^ )(A))/T 



-0 (1) (A) 



j^ r) f(X-X')Tln(l + e 



(0M(A)+r,M(A))/T 



) 



(r) dX'f(X-X') V ^(X') + Ef\X). 



\\\>x% 

(V) 

has been made. Here E\ (A) stands for the integral 



E 



(r) 



(A):= / /(A-A')Tln(l + e 



|0M(A)|/T 



(35) 



(36) 



and / represents the subscripts h and g. For T — > 0, the major contribution to the 
integral is from the regions near the zero points of (f)^(X), i.e. A = ±Ag r \ Therefore, 
the expansion of <p^(X) around A = ±Ag is 

W (A)=t« (A-Ar)+0[(A-Ar) 21 



where 



r) 



dX 



A=A, 



(r) 



(37) 
(38) 



Then the leading term of Ej (A) becomes 



Ef\x) 



2T 2 

7l 2 T 2 

6tM 



f (A - A?*) + / (A + A^ r) )] r rfx In (l + e -) 



'/ (A - 



(39) 



/ (A + A? )X 

Substituting these results into (1341) . we obtain the equations for <p r ) and r]( r \ namely 
(1) (A) = H- 2nP Sl (X) + h * (1)+ (A) + g * (2)+ (A), (40) 
(2) (A) = H - 2nPs 2 (X) + h * (2)+ (A) + g * (1)+ (A), (41) 



r? (1) (A) 



|A|>A« 



+ 
+ 



7l 2 T 



7T 



2rp2 



h(X- X')r) (1 \X')dX' 

h(X- A^) + h (A + A^ 
,(A-A«)+,(A + A«) 



(i) 



h{X-X')r) {2 \X')dX' 



(42) 



^ (2) (A) 



+ 
+ 



|A|>A< 2) 
7l 2 T 2 



|A|>A< 2) 



6t( 2 ) 

7T 



fr(A- A')?7 (2) (A')dA' 

/i (A - A< 2) ) + (A + a[, 2) 

^(A-A^)+,(A + A^) 

When /J 0, A{, r) -»■ oo, Q3D) and (jUD can be simplified to 
(r) (A) = H -2txPs t {X). 



h(X-X')r] m {X')dX' 



2rp2 



6f« 



(43) 
(44) 
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Since (r) (A{, r) ) = 0, s r (X { r) ) can be written as 

*r(4 r) ) = ^p- (45) 
From (JT3J) and ([14} we find that in the limit X { r) > 1, 

s r {tf) = ^e-*W'* + o(±). (46) 
Hence from both equations ([475]) and (|46p we obtain 



(r) = _ 3c 
2tt 



2tt 

where m 1 is an integral constant. 

Instead of (r) (A) and ?7 (r) (A), it is easier to work with the functions 



(47) 



S (r) {x) = ( e 2 ^ r) /3c K M M (A + A M ); for A > 0; 

| 0, for A < 0. 

rW(A) = { V r) (A + A^), for A > 0; (4g) 

[ 0, for A < 0. 

Notice that (j)^(X) > for |A| > Aq . Instead, we shift the integration variables so 
that they run from to oo. Observe that h(X + 2A r ' ) ) and g(X + X^ + A^) vanish as 
H -»■ (Xt ] -»■ oo) for finite A. The functions g(X + X { q ] - A? } ) and #(A - A^ + X { 2) ) 
remain finite. Equations P0l) - fj43l) in the form of the newly introduced functions are 



5 (D (A) = ^£ fl - K We-^) + - AO-S^AOdA' 

V3c v •'0 

+ /°° (/(A + Aj, 1} - A' - A^S^AOdA', (50) 
Jo 

T<1> < A > - wm + g(A "^w A ° 2>) + f A(A - y)T(1)(A ' )dA ' 

+ /°°9(A + A™ — A' — A< 2 >)T< 2 >(A')<iA', (51) 
Jo 

5(2) (A) = ^ fl- ^ 2 V 27rA / 3c ) + f°/i(A - X')S^(X')dX' 
V3c v 7 Jo 

+ r^A + ^-A'-AS^^CAOdA', (52) 



+ r g{X + X ( Q 2) - X' - AS°)T«(A')dA'. (53) 
jo 

Writing the equations in standard Wiener-Hopf form for — oo < A < oo gives 
S (1) W = f ( g\X) + b ( s\\) + r MA - X')S^(X')dX' 



+ r ^(A + Af ) - A' - A^S^A')^', (54) 
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T«(A) = f£\X) + b { ? (A) + r h(X - A')T«(A')<2A' 



+ 



g(X + X^-X'-X^)T^(X')dX', 



S( 2 )(A) = ff ) {X) + b ( s\X) + I" h(X - X')S {2 \X')d\' 



+ 



g(X + Aq 2) - A' - X^)S ( - 1 \X')dX', 
T^(A) = f^(X) + ^ 2) (A) + [°° h(X - X')T^\X')dX' 

J — oo 



+ 



g (X + xW-X>-XP)TM(X')dX', 



where 



4 r) (A) 

4 1} (A) 
^(A) 
/? } (A) 

/f(A) 
b?(X) 
#(A) 



2nP 





(r)„-27rA/3c 



for A > 0; 
for A < 0. 



0, for A > 0; 

-h * 5«(A) - 5 * S^(X + X ( o } - A{, 2) ), for A < 0. 

0, for A > 0; 

-h * SM(X) - g * S«(A + X ( o ] - Aj, 2) ), for A < 0. 

+ 5( 2 )'(o) ' tor A > u > 

for A < 0. 







MA) , g(A+A< 1) -A< 2) ) 

5< 2 )' (0) ~r S(!)'(0) 

0, 

0. 



for A > 0; 
for A < 0. 



for A > 0; 
for A < 0. 

for A > 0; 
for A < 0. 



(55) 
(56) 

(57) 

(58) 
(59) 
(60) 
(61) 
(62) 

(63) 
(64) 



—h * T«(A) - g * T( 2 )(A + A^ - X^] 
0, 

-h * T( 2 >(A) - g * T«(A + aE, 1} - X^] 
The equations fl54l) - (!571) can be solved via Fourier transform. Define the Fourier 
coefficients of 5 (r) (A) and T^(A) as 

dXe iXui S {r \X), 



— oo 
oo 



S {r \tu) -- 

f ir \u) = I" dXe iXw T {r \X). 

J — oo 

The functions S^(u) and T^(u) are analytic on the upper-half-plane. 
For simplicity, we write equations (!54|) - (l57|) in matrix form 

/oo 
K(A - A')S(A')dA', 
-oo 

/oo 
K(A - \')T(\')d\', 
-oo 

where 



S(A) 



/ SW(X) 
[ 5 (2) (A) 



T(A) 



/ t«(a) 

^ T&(X) 



(65) 
(66) 

(67) 
(68) 

(69) 
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b.(A)«(&, (70) 



#(A) 




f T(A) = | V(2> V; I > b T (A) = i i; 2) ^ i , (7i) 



k(a) = | m *(A+Af J -Any (72 ) 

1 ^(A + A( 2) -A«) MA) ; 1 J 

We denote Fourier transforms of S(A) and T(A) as 

S(u) = / d\e iXuJ S(\), (73) 

J — oo 

/oo 
dAe iAaJ T(A), (74) 
-oo 

which represents a Fourier transformation on each of the matrix components. 
Thus the Wiener-Hopf equations for S(A) and T(A) in Fourier space are 

S(u) = f s (u;) + bsM + £(w)S(w), (75) 
f(w) = f T (w)+b T (w)+K(w)T(w), (76) 

where 

Notice that when a; and A^ are real, the kernel K(w) is Hermitian, i.e. 

K\u) = K(u). (78) 
Equations f J75|) and fl76|) can also be written as 

S(co) = (I - K(w)) _1 (fsH + b s (w)) , (79) 
f (cu) = (i - K(w))" 1 (f T (w) + b T (w)) . (80) 

Since the function (l — K(w)) is nonsingular, hermitian and positive definite at 
u — 0, it is also positive definite for — oo < u < oo. From Theorem 8.2 of Gohberg and 
Krein [50], the function can be factorized as follows, 

(i-KHj^G+HG^M, (81) 

where each element of the 2x2 matrix G + (u) (G_(w)) is analytic on the upper-half 
(lower-half) plane. The limit G±(u) is chosen to satisfy the condition 

G±(oo) = I, (82) 
where I is the identity matrix. From (178)) . it follows that 

G+(-w) = G T M, (83) 
where A T represents the transpose of the matrix A. 
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In Appendix A and Appendix B, we have worked out the solutions to equations 
9]) and f lHUj) in detail. We showed in equations flA.9j) and (IB.9j) that their solutions are 



S(uj) = i^£ ( — 1 — 7 - ] G + (u;)G_(0) ( 1 ] (84) 

and 

TV) = (G + (w) - I)V (85) 
where V is a 2 x 1 matrix given in equation ( 1B.5|) . 



4. Universal low temperature properties and spin-charge separation 



To solve the TBA equations (J6J) and ffTOj) . we must find a way to express them involving 
the functions S and T. Let us rewrite the dressed energy from the first set of TBA 
equations © as 

Si(k - A) ln(l + e^i /T )ciA 

-oo 

s 2 (A;-A)ln(l + e^ /T )d\ 

-oo 

^k 2 -fi- 2ixP Sl * ai (k) -T I ^sAk-X) ln(l + e^ 1)/T )d\ 

J\x\>x ( l) 



- T f ° Sl (A; - A) ln(l + e^ )/T )dA - T / s 2 (£; - A) ln(l + e^ /T )d\ 
J—\ J\\\>\ 

A (2) 

-T ° s 2 (k-X)\n(l + e^^ T )d\ 
J-\ 

k 2 -fi- 2nP Sl * aAk) -T f sAk - A) ln(l + e^ /T )d\ 



- T U»« k - x ™ l+ ^^-£t%^ (86) 



where we evaluated the the integration between the Fermi points as 



2 T / ° Sr (A; - A) ln(l + e^ ' T )d\ « £ T /°° s r (A; - A) ln(l + e'^^A 

r=l "'-V r=l 

From the third line to the fourth line, we made use of the fact that 

SM'(O) = e^W^M. (88) 
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The other two integrals can be simplified as 



(r) 



r Sr (k - x - x^) 

Jo 



2vrA< ) r) /3c 2 T 2 (r) 

* (r) (A)+ Ll^ iX) 



K 



(r) 



+ 



oo e 27r(fe-A-2A< i r) )/3c 
47T 2 P 2 

(, loo 8tt 3 P 2 



S (r) (A)dA + 



6e -27rA^/3c 

00 7r 2 T 2 K (r:) e 27r(fc - A)/3c 



dA 



M*-A)/3c jS '(r)( A ) dA + 



6V3c 

oo 7r 2j.2 K (r) e 2 7 r(fc-A)/3c 



6^ 



T (r) (A)dA 

T (r) (A)rfA 





00 



e 2.(fc-A)/3 Ce -iuA£(r) (w)rfa;rfA 
2^(fc-A)/3c e -icAy(r)^^ a;rfA 



8vr 3 P 2 iu + 2vr/3c 



ttT 2 ^) T^H e27rfc/3( 
12^0 iw + 2vr/3c 



4vr 2p2 

Thus the dressed energy is 



y/3cP 2 /^) ~ (r) /2vri\ ttW ^ /27ri\ 



V 3c 



V 3c/ 



da; 



i9) 



2 -2 T 2 K (r) 



; /-oo ; IT I K" 

e(k)^k 2 -f,- 2nP Sl * ox(fc) - 2 £ jf Sr (fc - A) feM + ^W]dA - E ^7^7 
= k 2 — fi — 2-nPsi * ai(k) 



E 

r=l 



V3cW^)^ (r) /27TZ\ + 7T 2 T 2 K^ f(r) /2tTZ\ + 7V 2 T 2 K 



2vr 2 P 2 



\3cJ 3V3c V3c/ 3V3&SW(0)_ ' 



(0) 
(90) 



From (1A.9j) . S (27ri/3c) is expressed as 



(2^ = V3 ( 2_ni 



\ 3c 



V 3c 



G_(0) 



(91) 



To evaluate the sum Y%=i K ^ rS> S^(27ii/3c), we only need to consider the matrix product 

(«« &)&(™) = ^ P ( K n .( 2 )) G+ (f) G _(o)(;) 



= ^P(l 1 ) G + (0) G_ (0) ( J) 

= v^P (92) 

We made use of equations (1531) and (1A.8I) to derive the second line from the first line. 
The term G + (0)G_(0) was explicitly derived from its original definition ( 1HTT) . To do 
that, we needed to find K(0). Since 



*i(0) 



s l (x)dx = -, s 2 (0) 
o 



s 2 (x)dx = -. 

3 



(93) 
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hence 

v _i / 9 -1 \ 

(94) 



G + (0)G_(0) = (l-K(0))" 1 =(_ 2 i 



Similar to 



The next step is to evaluate the £r=i /t (r) [T^(2m/3c) + I / 3^(0) 
the previous calculation, it is easier to work with the matrix representations. Doing so 
yields 



(«a> ««))(t(^) + v) = («« ««) G+ (ff 

1 1 ) G + (0)V 



V 



4vr 2 P 
3\/3c 2 



S' r (0)V 



Here 



/oo 
si(/c — A)ai(A)dA 
-oo 



Si(X)ax(X)dX 



— oo 



1 In 3 



(95) 



2vr 2 P 

where we made use of the transpose of equation ( 1A.11I) . which is 

s,T(o)= S (i i)G+(o) - (96) 

Finally, the dressed energy can be written as 

3cP 2 cT 2 

£{k) nk 2 -ii- 2vrP Sl * ai (k) - — - — . (97) 



(9£ 



6^c 2ttc' 

Through integration by parts, equation ffl6|) is given by 

P=ir ^ £ 09) 

Trio l + e(^)/ T ' 1 J 

where the dressed energy is rewritten in the form 

e(k) = k 2 - A(T,H) (100) 

with 

3rP 2 rT 2 

A(T, P) = /i + 2vrP Sl * ai (k) + — + — . (101) 
For finite temperature, the pressure can be given by polylogarithm function 

F = -/5 riLi i ( - e * ) (102) 

This equation describes the exact low temperature thermodynamics and thus the full 
phase diagram can be given in finite temperature in a much larger regime than that 
given by field theory. However, the low temperature properties cannot be seen directly 
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from the polylog functions. Therefore, we apply Sommerfeld's expansion technique for 
this equation and study the universal low temperature properties. The expansion of the 
pressure is in terms of powers of (T/A(T, H)). Consider the leading order, the equation 



is 



2 A 3 - 

— 

3tt 



1 + 



7T- 



T 



8 \A 



(103) 



Furthermore, with the relation n = dP/dfi, after the iteration, and neglecting higher 
order terms of T 2 , H 2 and 1/c, the chemical potential fi is obtained as 



, oT l&Trnsi * aAk) 9cH 2 , 

„ = „V[i ^ + — (1 

3cT 2 T 2 
[1 + 37rnsi * ai(k)) + 



37TOSi * di(fc)) 



2n% 4 v v " \2nH 2 

Substituting ( 1104(1 into ( 11031) . the pressure P becomes 
2 



(104) 



P 



-n 3 ir 3 



8\cH 2 

1 — 6irnsi * ai(k) H t- (1 + 47TOSi * ai(fc)) 



571 7T 



3tt 

27cT 2 T 2 
+ „ c ; (1 + Aimsi * ai(k)) + (1 + 2imsi * ai(k)) 



(105) 



8n 5 7r 4 v v " 4n% 2 

Finally, with the relation F = fin — P, the free energy suggests a universal low 
temperature behavior of TLL, namely 

F = fin-P 



1 9cH 2 

-n 3 7r 2 (1 — ATlnsi * a>i(k)) „— 7 (1 + Qixns\ * aUk)) 

3 4n 2 7r 4 

3cT 2 T 2 

[1 + 6nnsi * Oi(fc)) (1 + 47rnsi * ai(fc)) 



4n 2 7r 2 

„ vrT 2 fC s C c 
Eq- — (— + — 

6 V v s v r . 



12n 



(106) 



where 



E 



1 9cH 2 
-n 3 7r 2 (1 — Airnsi * ai(k)) — - 9 _ d (1 + 6nnsi * ax(k)) 
3 



1 o / 27m 
—nix 1 - 



4n% 4 
2nln3\ 9c# 2 



3v/3c 



4n 2 7T 4 



1 



7m 



3nln3 



4 

— n 2 7r 3 (1 — 67msi * aAk)) 
9c 

4 9 o / 7m 3nln3 x 

— n 2 7r 1 = 

9c V v^c c / 

2n7r (1 — Annsi * ai(k)) 

( 2nn 2nln3\ 
2wn 1 



(107) 



(108) 



(109) 



\ 3V3c c J ' 
The spin and charge velocities can be derived from the relations v c = 
e' (k ) /2rcp c (k ) and v s = (j)^' (A ) /2rrp s (A ) [48] . For three-component fermions, there 
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Figure 1. (Color online) The left plot shows Cy — T curves for external magnetic field 
H = 0.001, chemical potential fj, = 6 and coupling constant c = 30. The right plot 
shows S — T curves for external magnetic field H = 0.001, chemical potential \i = 3 
and coupling constant c = 50. See text 

are two spin velocities v s i and v S 2, where v s \ = v S 2- The central charge for the spin part 
is C s = 2 and that for the charge part is C c = 1. The reason for C s = 2 is because 
the SU(3) invariant fermion model has two spin "Fermi seas" whose dependence on H 
are equal i.e., we considered the case where H\ — H 2 — H . This result shows that 
spin-charge separation exists for low temperatures and a small external magnetic field. 
This result is coincident with the results from conformal field theory [T2| [TT] . 

To compare the results from Sommerfeld's expansion and the polylog function, 
we plot graphs of specific heat Cy — T and entropy Sy — T versus temperature in 
FIG. [U Both graphs show that at low temperatures in the strong coupling regime, 
Sommerfeld's expansion agrees well with the polylog function. However, when the 
temperature increases, the deviation between the two curves become more apparent. 
This means that the Tomonaga-Luttinger liquid form of the free energy is only valid 
at very low temperatures and conformal invariance breaks down as the excitations take 
place at higher temperatures. 

5. Conclusion 

We have presented a systematic way to study low temperature behavior of systems 
with high spin symmetry and with repulsive interactions via the Wiener-Hopf method. 
In particular, we have derived the universal thermodynamics of ID strongly repulsive 
fermions with SU(3) symmetry under a weak external magnetic field. We have applied 
the Wiener-Hopf method to obtain the universal low energy physics in terms of spin- 
charge separation. The chemical potential (I104p . pressure f 1 1 5 j) and free energy f |106|) 
of the system have been derived analytically at low temperatures. The free energy 
gives insight into understanding the spin-charge separation theory for the higher spin 
systems. Moreover, beyond the spin-charge separation regime, the thermodynamics 
( I102j) has been given in terms of the polylog function. This result can be used to test 
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universal scaling behavior in the quantum critical regime. Our results open up further 
study of quantum criticality in ID interacting Fermi gases with higher spin symmetry 
and repulsive interaction. 
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Appendix A. Solution to Equation ( I79j) 



With the product decomposition given in (18111 . equation ( 1791) can be written as 

G;V)SM = G_(w)f s (w) + G_(w)b s (w). (A.l) 

The function G^ 1 (w)S(a;) is analytic on the upper-half plane by the definitions of 
G^ 1 (a;) and S(uj). On the other hand, G_(u;)bs(w) is analytic on the lower-half plane. 
We are left with the function G-(u)fs(w) which we have to decompose into the sum 
of a function which is analytic on the upper-half plane and another function which is 
analytic on the lower-half plane. 

First, notice that f$\uj) can be written as 



('■), 



2vrPi 



K 



(r) 



UJ 



a/3c \w + ie uj + 2m /3c J 
where e is small. The function G_(w)fs(w) in matrix form is then decomposed into 



(A.2) 



G_(w)! s (w) 



2vrPi 
VSc 



_i_( G _M-G_(-i £ ))(; 



(G_(w)-G_(-27ii/3c)) 



2vri/3c 



+ 



2vrPi 
73c 



1 



uj + le 



-G r-ie) 



27ri/3c 



G_(-2vri/3c) 



Thus ( 1A.1I) can be rewritten as 

G; 1 (a;)S(a;) - $+(u;) = + G_(w)b s (w) (A.4) 

where the left-hand side is analytic on the upper-half plane, and right-hand side is 
analytic on the lower-half plane. When uj — > oo, both the left-hand side and right-hand 
side tend to zero, thus 

G; 1 (w)S(w)-$ + (u;) = (A.5) 

which gives 

2vrPi 



SM 



G+(w)$ + (w) 

1 



UJ 



27ri/3c 



V3c 

G+(o;)G_(-27ri/3c 



G+ (o;)G_(0) 

uj + le 



K 



(1) 
(2) 



(A.6) 



Since 5^(0) = and S^ r \oj) is analytic on the upper-half plane, the following 



result 



lim u;S(u;) 

\uj\ — >oo 



is given by contour integration. This boundary condition is equivalent to 



G_(-27ri/3c) 



.(2) 



G_(0) 



(A.7) 



(Ai 
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Using this result we can show that 

2nP / 1 

V^c V^ + ie lo + 2m/3cJ ~ + v ~ 7 v " 7 ^ 1 



S(u;) = i^ (-^— - - , ^jjA G + (a;)G_(0) (]). (A.9) 



Besides, note that 



S'(0) = ^-S(A)| = - lira c 2 S(u;), (A.10) 

thus we find the expression of S'(0) 

4vr 2 P / 1 \ 

S(0 »^ G - (0 >(i)' (A ' n) 

Appendix B. Solution to Equation (1801) 

Similarly for T(ui) with the factorization (IHTj) . we find that 

G-\u)T (w) = G_(w)f T (w) + G_(w)br(w). (B.l) 

With similar arguments as before, we only need to find a decomposition of G_(w)fr(a;) 
into the sum of two functions that are analytic on the upper-half plane and lower-half 
plane, respectively. Consider the functions 

f (D m fe (A) | ^A + A^-A^) 
f(2)m MA) g(A + A^-Ag) 

/t (Aj ~ ^(0) + Who) • [ti - 6) 

After taking a Fourier transformation, the matrix fr(w) can be expressed in terms 
of the kernel K(w) as 

f T (w) = K(w)V = (i - Gl 1 (o;)G; 1 (a;)) V, (B.4) 

where 

V - f 1/5(1) ' (0) 1 (B 5) 

Multiplying G-(u>) on both sides of this equation, we have 

G_(w)f T (w) = (G_(w) - G; 1 ^)) V. (B.6) 
Substituting this relation back into equation (IB. ip gives 

G; 1 (o;)f (w) + G; 1 (w)V = G_(w)V + G_(w)b T (w). (B.7) 

Based on analyticity arguments, matrix elements on both sides of the equation are 
equal to a function that we shall denote in matrix form by Q(w), i.e., 

G; 1 (w)T(u;) + G; 1 (u;)V = QH. (B.8) 

Similarly to S(u), the function T(w) vanishes when u; — > oo. Because Q(w) is an entire 
function that is bounded, it is equal to a constant by Louville's Theorem. This constant 
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can be determined by taking the limit u — > oo which gives Q = V. Hence we finally 
obtain an equation for T(u) expressed in terms of functions whose values are known, 
i.e., 

f(u) = (G+(w) - I)V. (B.9) 
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